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ABSTRACT 

In recent years there has been much interest in the nature of two stars, V407 Vul 
and RX J0806+1527, which are widely thought to be binary white dwarfs of very 
short orbital period, 570 and 321 seconds respectively. As such they should be strong 
sources of gravitational waves and possible ancestors of the accreting AM CVn stars. 
Monitoring at X-ray and optical wavelengths has established that the period of each 
star is decreasing, at rates compatible with that expected from gravitational radiation. 
This has been taken to support the "unipolar inductor" model in which the white 
dwarfs are detached and the X-rays produced by the dissipation of magnetically- 
induced electric currents. In this paper we show that this interpretation is incorrect 
because it ignores associated torques which transfer angular momentum between the 
spin of the magnetic white dwarf and the orbit. We show that this torque is ~ 10 5 
times larger than the GR term in the case of V407 Vul and ~ 10 times larger for 
RX J0806+1527. For V407 Vul, the unipolar inductor model can only survive if the 
white dwarf spins ~ 100 times faster than the orbit. Since this could only come 
about through accretion, the validity of the unipolar inductor appears questionable 
for this star. We also consider whether accretion models can fit the observed spin- 
up, concluding that they can, provided that a mechanism exists for driving the mass 
transfer rate away from its equilibrium value. 
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1 INTRODUCTION 

Over the past decade, observations have established the 
existence of a population of some 10 0-200 million dou- 
ble white dwarfs with in our Galaxy jMarsh et alJ Il995t 
iNapiwotzki et alJ 12003^ . These and their descendants are 
thought likely to be a dominant source of low frequency 
gravitational waves in the Galaxy (e.g. iHils et all Il99dt 
iNelemans et al.ll2001bf) . and are a possible progenitor popu- 
lation of Type la supernovae. A significant fraction of these 
binary stars are close enough that gravitational wave losses 
will cause them to undergo mass transfer within a Hub- 
ble time. Most will merge to form single stars, variously 
suggested to be Type la supernovae, or in the major i ty of 
cases , sdB, sdO or R CrB stars (e.g. IWebbinklll98i liber] 
Il990l : ISaio fc Jeff ervl 1200^ . If any systems survive the on- 
set of mass transfer as binary stars, then they would be- 
come semi-detached accreting binary stars with white dwarf 
donors, properties matched by the AM CVn stars, which 
feature helium-dominated spectra and orbital periods which 
range from 10 to 65 minutes. 

The question of survival as a binary is key to whether 



double white dwarfs are the ancestors of AM CVn stars and 
is important to the prediction of gravitational waves. The 
nearest we can get to a proof that the start of mass trans- 
fer can be survived is to identify accreting pairs of double 
white dwarfs with periods short of 10 minutes since there 
are a lternative routes for sy s tems with periods longer t han 
this llNelemans et alJl2001at iPodsiadlowski et alJl200&t . In 
this context, two stars, V407 Vul and RX J0806+1527 
have generated much interest in the past few years because 
they show periods of 570 and 321 seconds respectively, and 
there are several rea sons to think that these may be orbital 
llCropper et al Il99gft . rather than, for example, the spin of a 
magnetic accr eting white dwarf a s seen in "intermediate po- 
lars" (but see lNorton et al.ll2004T) . The key features of these 
systems are (i) X-ray light curves which are off for about 
half the period llCropper et aljll99il : llsrael et al]ll999T) , sug- 
gestive of a spot of emission on a spinning star, (ii) opti- 
cal and i nfrared light curves that show the same perio d and 
no other iRamsav et all200fll2002al llsrael et all20oA. (in) 
weak optical line emission iRam^av^^dT2^021anisrael et alJ 
2002), and (iv) very soft X-ray spectra. Items ii-iv in partic- 
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ular count against TNorton et alJ i2004T) 's intermediate polar 
model. 

Even on the assumption that we are seeing orbital peri- 
ods, questions remain over the nature of these stars. There 
are currently three popular models which are (a) the po- 
lar model in which we see X-rays from the magnetic poles 
of a white dwarf locked to its companion dCropper et alJ 
1998), (b) the unipolar inductor (UI) model in which X- 
rays comes from the dissipation of electric currents gen- 
erated from a slight asynchronism betwe en the spin of a 
magnetic white dwarf and its companion dWu et alJl200l . 
and (c) the direct impact model in which X-rays come 
from the point at which the mass transfer stream impacts 
the white dwarf, which can happen directly in these very 
compact system s dNele mans et al.ll2001at iMarsh fc Steegha 
120021 ; iRamsav et alJl2002rJ) . The polar and direct impact 
models rely on accretion and so imply that the systems have 
survived contact; the system in the UI model is not in con- 
tact and therefore has not necessarily had to survive it. The 
polar model has fallen somewhat out of favour because of 
the weak optical line emission, an absence of or, at best, 
weak polarisation iRamsav et all2002bl : lReinsch et"al .12004 
llsrael et alJl2004) and the X-ray spectra, which are u nusu- 
ally soft, even for polars (but see lCropper et al]l2004|) . The 
two most promising models left are the UI and direct im- 
pact models. It was soon realised that period changes might 
help distinguish between these models since on the detached 
UI model it was predicted that the period would decrease 
while in accreting models the period should increase as the 
donor expands. In a series of papers, the periods of both 
V407 Vul and RX J0806+ 1527 have been definitively es- 
tablished to be decreasin g dStrohmaverll2002l : Hakal a et al 
| 2003|: IStrohmaveJ pool llsrael et alJ 12004 Hakala et a" 
l2004l : lRamsav et al.l200alStrohmaverl2005|). which, in spite 
of some problems which f ace the UI model JlVlarsh fc Steeghsl 
120021: iBarros et alJl2005^. has been taken to be strong evi- 
dence in its favour dCropper et"alll2004) . 

In much of this work it has been assumed, almost by 
default, that, since the systems in the unipolar inductor 
(UI) model are detached, their orbital periods change at 
a rate governed by gravitational radiation alone (see e.. 
StrohmaveiH2002l: Icropper et aljEool | iakal a et al 1200 



Willes fc Wull200^ . In this paper we show that this is very 



far from being the case, and th at on the contra ry the ob- 
served period changes rule out IWu et all |2002) 's original 
model once- and- for- all. We further show that it can only 
survive in a very different form which necessarily requires 
accretion to have occurred, raising the question of why it 
ever should have ceased. We begin by discussing the period 
change expected for the UI model. 



2 PERIOD CHANGES FROM INDUCTION 
TORQUES 

The key point about the UI model which has not always 
been appreciated is that although the two white dwarfs are 
detached, their period evolution reflects the combination of 
angular momentum loss through gravitational radiation and 
induction-driven angular momentum interchange within the 
system . Although this was stated explicitly by IWu et al] 
(2002), who also present equations which correctly account 



for the effect, its importance seem s not to have be en gener- 
ally recognised. As we will show, in lWu et alJi2002f) 's model, 
the angular momentum interchange term dwarfs the GR loss 
term, and therefore the idea that the UI model evolves in the 
same way as a pair of detached stars is completely wrong. 

In the UI model the X-rays come from the dissipation 
of electric currents generated by magnetic induction. This 
implies that a torque T exists which transfers angular mo- 
mentum between the spin of the magnetic primary star and 
the binary orbit. If the spin angular frequency of the pri- 
mary star is fl„, then the rate of work done by the torque 
is TQ S (we define the sign of T such that a positive value 
spins the primary star up). An equal but opposite torque 
acts to extract angular momentum from the orbit, removing 
energy from it at rate TQ where f2 is the orbital angular 
frequency. The difference between the two is dissipated, at 
least in part, in the form of X-rays so that 



L x <T(n o -fl s ). 



(1) 



Therefore the rate of change of orbital angular momentum 
J is given by 

Lx 



J = Jgr 



(2) 



(l-a)JV 

where we assume from now on that equation is an equal- 
ity and where we have introduced a = Q 3 /Q , following 
IWu et al.N2002l) . In this equation Jgr is the rate of angular 
momentum lost due to gravitational waves alone. For two 
detached stars, the rate of angular momentum loss imme- 
diately leads to the period change through P/P — 3 J/ J, 
therefore whether the period change occurs at the GR rate 
depends upon the relative magnitudes of the two terms in 
Eq. [5] If we multiply through by the orbital angular fre- 
quency, then we obtain an illuminating version: 

/X (3) 



E = —L 



GR 



(1-a) 



where Lgr is the luminosity in gravitational waves and E is 
the rate of change of orbital energ y. Equation El is i dentical 
in physical content to Equation 9 of lWu et al.l d2002l) . except 
that for simplicity we have ignored the moment of inertia of 
the secondary star as it makes only a small difference to the 
results. Equation |3| shows that the simple GR formula for 
the period derivative can only be used if 

L x <|1- a\L G R- (4) 

We now show that, on lWu et alJ (|2002)'s model, this is not 
at all the case for either V407 Vul or RX J0806+1527. 



2.1 V407 Vul and RX J0806+1527 

In order to estimate the gravitational wave luminosities, let's 
assume that we are seeing the pure GR rate, in which case 
it can be shown that 



Lgr = 



5c 5 P 2 



(5) 



1152Gtt 2 ' 

independent of the masses. The observed rates of pe 



3 x 10~ 1J s 
3.7 x 10" 11 s/s 



the two systems. 
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The X-ray lu minosity of V407 Vul is ~ 
lO^Cri/lkpc) 2 ergs/s (iRamsay et al.ll2005l) . From Eq. 
and assuming IWu et alJ (|2002h 's value of 1 - a = 0.001, 
one would therefore predict that the induction torque term 
is ~ 10 5 (d/l kpc) 2 times larger than the GR term. In other 
words the observed rate of period change in V407 Vul is far 
too sm all compared to a prediction based upon IWu et all 
(2002)'s model. The model is nearer the mark in the case of 
RX J0806+1527 fo r which L x ~ 2 x 10 33 (d/l kpc) 2 ergs/s 
iStrohmaved l2005t) . and the predicted induction torque 
term is "only" ~ 10 times the GR term. However, this still 
shows that it is incorrect to assume that period changes in 
the UI model run at a rate given solely by GR. 

A problem for the UI model has always been its short 
lifetime befor e synchron i sation occurs, estimated to be ~ 
1000 years bv lWu et alJ (120021) . In order to have (at least) 
two such systems within a kiloparsec of the Sun, this sug- 
gests a formation rate of order 1 per year within the Galaxy, 
unless the asynchronism can be regenerated by some as- 
yet-undetermined mechanism dCropper et al .1200411 . We will 
now show that the problem is actually much worse than even 
this estimate suggests. During synchronisation, the spin rate 
of the primary will change much faster than that of the orbit, 
so we can estimate the timescale by neglecting the change of 
orbital period. The rate of change of the spin energy is then 
given by the power injected into the spin = TO s so using 
Eq.HJ 



1- £ Q 



Q s Lx 



(6) 



where I ~ (1/5) Mi R\ is the moment-of-inertia of the pri- 
mary star. Dividing through by f2 Q gives 



1 



IQ% 1- a' 



(7) 



This equation is identical to equation 8 of IWu et alJ ( 200 
who also find that Lx oc (1 — a) 2 , so we can write 



T i L x) f .2 
L X = — —{t — Ct) 



(l-a)8 

where the subscript zeroes indicate initial values. Therefore 



mi{i-a)i {1 ~ a) ' 

which can be integrated to give 



(1 - q) e 



-t/T 



with the synchronisation timescale r s given by 
t s = — (1-a) . 



(9) 



(10) 



(11) 



Since Lx oc (1 — a) 2 , the X-ray luminosity decays with a 
time constant of r s /2. Taking Mi = 0.6 M©, Ri = 0.01 R©, 
we find r s ~ 4400(l-a) 2 (d/1 kpc)" 2 years for V407 Vul, and 
7 x l0 5 (l-q) 2 (d/l k pc)~ 2 years for RX J0806+1527. Putting 
in lWu et all (I2002T) 's value of 1 - a = 0.001 leads to X-ray 
flux decay constants of 0.8 days for V407 Vul and 4.2 months 
for RX J0806+1527 which are incompatible with the 15-odd 
years that these systems have been followed, let alone plau- 
sible formation rates. Our timescale for V407 Vul is a great 
deal s horter than the ~ 1000 years quoted by | Wu et alJ 
(2002). It is not possible to determine from IWu"et^dT "the 




log t (yr) 



Figure 1. Synchronisation in the unipolar inductor model. From 
top to bottom the panels show the degree of asynchronism, the 
period derivative scaled by the GR-only value and power dis- 
sipated as a function of time. The left-hand panel starts with 
1 - a = 0.001 and L x = 10 35 erg/s to match V407 Vul while the 
right-hand panel similarly matches RX J0806+1527. Note that 
the horizontal axes of the two panels differ from each other. 



exact cause of this difference, but we think that it is sim- 
ply because their estimate was not specifically for V407 Vul, 
but applies to a lower luminosity model. Whatever the rea- 
son, the important point is that the timescale is in fact 
much shorter than their 1000 year value, which was already 
hard to credit. In Fig. we show numerical integrations of 
the synch r onisat ion for fixed dipole moments according to 
IWu et all (l2002T) 's equations. The initial conditions of the 
left- and right-hand panels were chosen to match V407 Vul 
and RX J0806+1527 respectively. The integrations confirm 
that the UI model can only sustain the observed X-ray fluxes 
for a very short time indeed. In the integrations we account 
properly for the changing orbital period, and hence in the 
right-hand panel a small residual asynchronism is seen at 
long times as the spin frequency must always lag the orbital 
frequency by a small amount. 



2.2 Revision of the Unipolar Inductor model 

We have shown that lWu et alJ ll2002l) 's model is ruled out by 
the measured rates of period change in V407 Vul and by its 
very short lifetime for both V407 Vul and RX J0806+1527. 
We now consider whether it can be adjusted to match them. 
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The problem stems from the inefficient X-ray generation as 
only 1 part in 1000 of the power transferred from the orbit 
to the spin is dissipated. This can be avoided if very different 
values of a = Q s /Q are considered. For V407 Vul, a value 
of |1 — a\ > 100 would reduce the induction torque term to a 
value comparable to or smaller than the gravitational wave 
term. Since we expect prograde rotation (a > 0), this means 
that a > 100, i.e. that the magnetic white dwarf must spin 
of order 100 times faster than the orbit so that P sp i n < 6 sec, 
which is or order the break-up spin rate of typical white 
dwarfs. Since a > 1, the induction torque acts in the reverse 
sense to GR and therefore the large value of a is needed not 
only to obtain the correct order-of-magnitude for P, but also 
to give the correct sign. 

Such a large value of a can only be attained through 
accretion-driven spin-up. The increased value of 1 — a also 
means that the magnetic mome nt of the white d warf must 
be much less than supposed bv IWu et al. ( 2002|) . Accord- 
ing to their model, the power dissipated in electric currents 
scales as oc ^i 2 (l — a) 2 , where fi is the magnetic moment of 
the white dwarf. Since in V407 Vu l we require 1 — a to be 10 5 
times larger than lWu et al ] d2002l) assume, the magnetic mo- 
ment of the white dwarf must correspondingly drop by 10 s , 
which would therefore end up at about 10 2r to 10 28 Gem 3 , 
corresponding to a surface field of only 1 to 10 G. We would 
then have a case of accretion onto a very weakly magnetic 
white dwarf (non-magnetic as far as one could tell from ob- 
servations), which begs the question of why the accretion 
should ever have stopped to allow the unipolar inductor to 
get going in V407 Vul but not in similar systems, such as the 
AM CVn stars and cataclysmic variable stars. Occam's razor 
suggests that if accretion has to be invoked to allow the UI 
model to work at all, then one should favour accretion-only 
models, which we will look at in the next section. 

The much lower X-ray luminosity of RX J0806+1527 
leads to less stringent requirements. For 1 — a = 0.001 the 
induction term is only 10 times the GR term, and so a mod- 
est increase leads to an acceptable match given the uncer- 
tainties. A tougher constraint comes from the X-ray decay 
timescale, which at 4 months for 1 — a = 0.001 is much too 
short. However, an increase to 1 — a = 0.01, gives a time 
scale of ~ 30 years which is probably long enough that it 
could not be ruled out by observations. 

The larger values of 1 — a needed to match V407 Vul 
and RX J0806+1527 have one further consequence: they 
make it easier to see the phase shifts between optical and X- 
ray pulses predicte d under the UI mod el, but not observed, 
which were used bv lBarro s et alj (l2005h to rule out the dipo- 
lar field geometry used bv lWrH^talj i2002h . For small values 
of 1 — a, it was always possible that an unlucky distribution 
of observations had lead to our missin g these phase shifts. 
The raised values oil — a and Fig. 6 of lBarros et al] (2005) 
show that this is no longer the case. 

In all the above, the X-ray luminosity, and therefore 
the distance, is crucial. Smaller distances reduce Lx and 
therefore alleviate the problems discussed above. We have 
assumed d — 1 kpc in each case. For V407 Vul this comes 
from assuming that the variable is at the same distance as 
the G star which dominates its spectrum (Steeghs et al., 
in prep.). For RX J0806+1527 on the other hand, the esti- 
mate comes from its blue colour, magnitude and adopting 
My = 11 as the absolute magnitude of a white dwarf of 



comparable colour. Neither constraint is secure, and lower 
distances are certainly possible. A crude lower limit comes 
from the absence of detectable proper motion in either star, 
which suggests that d > 100 pc. However, it is hard to see 
that RX J0806+1527 can be near this limit at the same 
time as being blu e and faint (U — B — —1.1, V = 21.1 
llsrael et alJ 12002ft . unless the emitting area is much less 
than that of an average white dwarf. A similar argument 
applies to V407 Vul, because although it is brighter than 
RX J0806+1527, much of its flux in V comes from the G 
star, leaving the variable comparable to RX J0806+1527 
in brightness. Thus d > 300 pc, is probably a better lower 
limit, so the X-ray fluxes could perhaps be lowered by a 
fact or of 10. This wo uld leave RX J0806+1527 compatible 
with lWu et"al"l(l2002r) 's model, but V407 Vul still a good way 
from it. Better constraints on the distances to these systems 
are needed. 

We now look at the whether accreting models can do 
any better in matching the observed period decreases in 
RX J0806+1527 and V407 Vul. 



3 PERIOD CHANGES IN ACCRETING 
BINARY STARS 

Assuming conservative mass transfer one can show that 

?-(7-<->£)- 

where q = M2/M1 and M2 is the mass of the donor. Since 
for stability we require that q < 1, and since J < and 
M2 < 0, mass transfer offsets the angular momentum loss 
term. The angular momentum term depends upon whether 
the accretor's spin is strongly coupled to the orbit or not. 
If it is we have si mply J = Jgr', if not then Eq. 1 from 
iMarsh et ail (I2004T) gives 

j = j G R + VGMiR h M 2 , (13) 

where Rh is the circularisation radius, which effectively in- 
creases the rate of loss of orbital angular momentum. 

An illuminating comparison with the section [5] can be 
made by calculating the "effective" angular momentum loss 
rate J e // such that the detached formula can still be applied 
as P/P = 3J e ff/J because this can then be compared di- 
rectly to Eq. |2] Assuming that L x = -GM1M2/R1 where 
Mi and Ri are the mass and radius of the accretor, and 
using Eq. 1121 f strong coupling), one can show that 

Jeff = Jgr + (l-q) — ^, (14) 
a I to 

in the strongly coupled case where a is the orbital sep- 
aration. Comparing with Eq. 121 we see that th e factor 
1/(1 - a) ~ 1000, which causes IWu et al] J2002ft 's model 
so much trouble, is replaced by (1 — q)R\/a ~ 0.1. This is 
a considerable improvement, and given the uncertainties in 
Lx and Lqr discussed in section |2~!H could be made consis- 
tent with the observed P values for both systems. 

For RX J0806+1527, the X-ray luminosity is already 
lo wer than the gravita tional wave luminosity, as pointed out 
bv lStrohmave T lElOSft . and so Eq. 1141 predicts that the ob- 
served period change must be close to the pure GR one. 
Weak coupling is better still since it increases the loss of 
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angular momentum from the orbit, but since it is otherwise 
qualitatively the same, we won't discuss it further. 

This discussion hides the ugly truth about accreting 
models which is that the equilibrium value of M2 leads to 
an increasing period. Put differently, the current mass trans- 
fer rates in V407 Vul and RX J0806+1527 must be < 60% 
of their equi librium values in o rder for their periods to de- 
crease at all jMarsh et all2004h (with some uncertainty over 
V407 Vul since some of its apparent period change could be 
caused by light-travel time effects if it is truly associated 
with the G star mentioned earlier). This is the price one 
must pay to accept accreting models. We regard it as a small 
one because there are already examples of systems which do 
not have mass transfer rates that match expectations, such 
as long period dwarf novae, and there are ways to make the 
mass transfer rate deviate from its equilibrium value as we 
will discuss below. 

If the mass transfer rate is below its equilibrium value, 
how long will it stay like this and would we expect to see 
a noticeable increase over the course of a few years? The 
mass transfer rate in the two systems is expected to b e in 
the adiabatic regime <Webbinkjl984 iMarsh et alJl2004h for 
which M oc A 3 where A = R2 — Rl is the amount b y 
which the donor overfills its Roche lobe jMarsh et alj|2 004). 
Therefore the rate of change of mass transfer rate is given 
by 

M ■ 

M = 3— A, (15) 

and so the timescale for significant change of the mass trans- 
fer rate is given by 



(16) 



_ M _ 1 A 
Tc ~ m" ~ 3 A' 

Using Equation 19 of jMarsh et al.ll2004T) . assuming strong 
coupling and neglecting the M dependent term in order to 
calculate the maximum rate of change of M gives 



Tc ~ = - 



1A_ 

&R~ 2 



J o r 



-Jgr 



(17) 



Taking Mi = 0.6 Mq, and M 2 = 0.07 and 0.12 M for 
V407 Vul and RX J0806+1527 resp ectively, and using equa- 
tions 10-12 of iMarsh et alJl2004) to calculate A, we find 
timescales of ~ 1000 and 100 years for significant alter- 
ations in the mass transfer rate to occur in V407 Vul 
and RX J0806+1527 respectively. These numbers suggest 
that below-equilibrium transfer rates can be sustained for 
many years, as observed. We must appeal to some un- 
known mechanism to force the departure from equilib- 
rium in the first place, but there is no shortage of can- 
didates. For example, a star-spot moving over the inner 
Lagrangian point is o ne possib i lity IILivio fc Pr inglc 1994; 
iHessman et alJ |2^0fjJ^l^s^mari[j200^), irradiation-induced 
cycles another lRltteret^dTl2000l) . widening of the orbit 
caused by the m ass ejected in nova explosions a third 
iShara e t aljl986 ]) and synchro nisation-induced detachment 
a fourth llLamb fc MeTIalll987ft . One can even circumvent 
the relation between mass transfer rate and period deriva- 
tive a ltogether through alterations in the structures of the 
stars ( Applegatel Il992l) . a mechanism which is thought to 
be responsible for cyclical variat ions in the periods of many 
eclipsing cataclysmic variables feaptista et al.l Eo03). This 



mechanism can make the period decrease even if mass trans- 
fer occurs at the equilibrium rate, although the low X-ray 
luminosity of RX J0806+1527 means that there remains a 
need for below-equilibrium mass transfer whatever the cause 
of the period decrease. Which of these mechanisms, if any, 
can be applied to double white dwarfs is not clear, but nei- 
ther is it clear that any of them are not applicable. 



4 CONCLUSIONS 

We have shown that, contrar y to widespread assumption, 
the unipolar inductor model llWu et all I2002T) does a poor 
job at predicting the magnitude of the period changes ob- 
served in the two candidate ultra-compact binary stars 
V407 Vul and RX J0806+1527 in the sense that the pre- 
dicted changes are much larger than those observed. This 
removes the main piece of evidence supporting the model. 
The reason for this is that if there is only a small asyn- 
chronism between the spin and orbital periods, then much 
more energy is transferred between the orbit and spin than 
is dissipated in X-rays. The problem is particularly acute in 
the case of V407 Vul for which the unipolar inductor model 
predicts a rate of period change some 100,000 times greater 
than is observed. The unipolar inductor model can only ap- 
ply to V407 Vul if its primary star rotates much faster than 
the 570 second putative orbital period; a larger asynchro- 
nism is also required to allow the unipolar inductor model 
to last more than a few years. This suggests that for the 
unipolar inductor model to work, accretion is a necessary 
precursor, without there being any obvious reason for it to 
have ceased. The problem with accretion-only models re- 
mains that the equilibrium mass transfer rates should lead 
to increasing orbital periods (as opposed to the observed 
decreases), and so the two systems must currently be trans- 
ferring mass at below equilibrium rates or other mechanisms 
must be affecting the orbital periods if accretion is to work. 
Given that there are many examples of systems where this 
is the case, this difficulty seems less significant than those 
facing the unipolar inductor model. Nevertheless, the nature 
of these systems remains far from clear, and observational 
efforts to pin them down must continue to be pursued. 
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